Observation of an in-plane vortex lattice transition in the multiband superconductor 
2i/-NbSe 2 using scanning tunneling spectroscopy 



m 
o 

(N 



I. Fridman, 1 C. Kloc, 2 C. Pctrovic, 3 ' 4 and J. Y.T. Wei 1 ' 4 

1 Department of Physics, University of Toronto, 60 St. George St., Toronto ON M5S1A7 Canada 

2 School of Materials Science and Engineering, Nanyang Technological University 639798 Singapore 

3 Condensed Matter Physics and Materials Science Department, 
Brookhaven National Laboratory, Upton, NY 11973 USA 
4 Canadian Institute for Advanced Research, Toronto, ON, M5G1Z8 Canada 

We report real-space evidence for an in-plane vortex lattice (VL) transition in the multiband su- 
perconductor 2//-NbSe2, observed by laterally imaging subsurface vortices with scanning tunneling 
spectroscopy. With magnetic field applied || afc-surface up to 6.5 T at 300 mK, spatial maps of the 
zero-bias conductance revealed a stripe pattern whose dominant periodicity showed a discrete shift 
near 0.7 T. The stripes can be interpreted as the surface projection of a distorted hexagonal VL that 
undergoes a first-order reorientation transition. This transition occurs for field || [100] but not [110], 
and is correlated with the multigap characteristics seen in the heat capacity and thermal transport 
data. These results implicate the multiband pairing of 2ff-NbSe2 in the VL transition we observed. 
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PACS numbers: 74.55.+V, 74.25.Uv, 74.70.Ad 

In the mixed state of type-II superconductors, a dia- 
magnetic current flows along the sample edge and para- 
magnetic currents loop through an ordered array of flux 
quanta in the bulk to form a vortex lattice (VL) [1]. The 
superconducting VL is known to be sensitive to crys- 
tal anisotropy, vortex- vortex interaction, pairing symme- 
try and band structure [2—5]. These factors can cause 
the VL to undergo geometric variations and thermody- 
namic transitions, which have been observed in various 
superconductors by both reciprocal-space and real-space 
measurements [6-15]. A transition of particular interest 
involves the rotation of a hexagonal VL vs. c-axis mag- 
netic field, observed by neutron scattering in the multi- 
band superconductor MgB2 [16]. This VL transition is 
second-order in character, and is correlated with the field- 
induced suppression of pairing in one the bands. In gen- 
eral, field depairing in a bulk superconductor is a function 
of its energy gap and Fermi velocity. When the pairing 
involves multiple bands [17-19] of different dimension- 
alities, as in the case of MgB2, the depairing can thus 
proceed band by band and depend on the field direction 
[20, 21]. Given the discreteness and anisotropy of this 
depairing process, it is natural to consider whether other 
types of VL transitions could occur in such multiband 
superconductors . 

In this Letter we report real-space evidence for a first- 
order VL transition in superconducting 2ii/-NbSc2, ob- 
served by scanning tunneling spectroscopy (STS) with 
magnetic field H applied parallel to the a&-planc. 2H- 
NbSe2 is a well-studied superconductor that has shown 
abundant evidence for multiband pairing [22-32] . Spatial 
maps of the zero-bias tunneling conductance over the ab- 
surface, measured at 300 mK and up to 6.5 T, revealed 
a field-dependent stripe pattern that can be identified 
as the lateral projection of a distorted hexagonal VL. 
While the stripe periodicity varies expectedly as 1/ VH, 



adjacent stripes show an amplitude alternation that in- 
creases with the field. Fourier analysis of the pattern re- 
vealed two peaks whose relative intensities are reversed 
near 0.7 T, corresponding to a discrete shift in the dom- 
inant stripe periodicity. This shift is observed for H || 
[100] but not [110], and is correlated with the multigap 
characteristics seen in field-dependent heat capacity [28] 
and thermal transport [22, 23]. We interpret this period- 
icity shift in terms of a VL reorientation transition that 
is driven by field depairing over parts of a multi-sheeted 
and anisotropic Fermi surface [26, 27, 33, 34]. These re- 
sults implicate the multiband pairing of 2_£Z-NbSc2 in the 
VL transition we observed. 

Our STS measurements were made on single crystals 
of 2_ff-NbSe2, using a horizontalized scanning tunneling 
microscope housed inside a 3 He refrigerator, with the 
field applied vertically by a superconducting solenoid. 
The crystals were grown by an iodine-vapor technique 
[35], structurally verified by x-ray diffraction (XRD), and 
showed superconducting critical temperatures sa 7.2 K. 
The crystals were ~ 5 x 5 x 0.5 mm 3 platelets with wide 
c-axis faces, and were aligned in-plane by XRD before 
being cleaved and cooled in zero field to 300 mK. Topo- 
graphic images showed atomically-resolved flat surfaces 
(see Fig. SI of the Supplemental Material [36]). The 
dl/dV conductance spectra were acquired using lock- in 
amplification with 20 fTV excitation at 505 Hz, and the 
typical junction impedance was 10 Mil at 4mV bias. Be- 
fore measurement, the Pt-Ir tips were field-emitted in 
situ to ensure junction stability over long scans, and RF- 
filtering was used throughout the cryogenic wiring to at- 
tain maximum spectral resolution. Reproducibility of our 
measurements was checked over different scan areas and 
on different crystals. 

Our cryomagnetic STS technique is illustrated Fig. 1. 
The STS geometry is shown in Fig. 1(a), with a diamag- 
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FIG. 1. (color online), (a) Scanning tunneling spectroscopy 
on the afe-surface of 2.ff-NbSe2 at 300 mK, with a field-induced 
diamagnetic screening current (red line) flowing across the 
surface and paramagnetic currents (blue loops) encircling sub- 
surface flux lines, (b) Spatial variation of the normalized 
dl /dV spectra, as the tip is scanned along a 140 nm line per- 
pendicular to the in- plane field at 0.1 T. A two-dimensional 
map of the normalized zero-bias conductance Go revealed a 
field-dependent stripe pattern as shown in Fig. 2a. (c) Profile 
plot of the Go map across the stripes at 0.45 T, after line- 
averaging along the field, showing the alternating amplitude 
between adjacent stripes. 

netic screening current sweeping across the sample sur- 
face, and paramagnetic vortex currents looping around 
flux lines in the bulk. In contrast to previous STS stud- 
ies which image cross sections of the VL [37-41], our 
STS technique images the subsurface VL laterally and 
with depth sensitivity, by virtue of the interaction be- 
tween surface and vortex currents [42]. Essentially, this 
interaction causes the field-dependent dl/dV spectra to 
vary spatially, as shown in Fig. 1(b) for data taken on the 
afr-surface of 2i/-NbSe2 along a line perpendicular to the 
in-plane field at 0.1 T. The spectra shown were normal- 
ized relative to the above-gap conductance at 3 mV. A 
planar map of the similarly normalized zero-bias conduc- 
tance Gg(r), shown in Fig. 2(a), revealed a stripe pattern 
that is quantitatively related to the VL spacing. Adja- 
cent stripes alternate in amplitude, as shown in Fig. 1(c) 
by the stripe-crossing profile of Go(r) at 0.45 T, plotted 
after averaging along the stripes. This amplitude alter- 
nation can be attributed to the depth alternation of the 
subsurface vortices, allowing us to associate the major 
(higher amplitude) and minor (lower amplitude) stripes 
with the first and second vortex rows, respectively. 

Identifying the stripe pattern as a lateral projection of 
the subsurface VL, we examine the field evolution of our 
Go data. Fig. 2(a) displays a selection of Go(r) maps 
for H || [100] in the range of 0.25 to 0.8 T, which is well 
into the mixed state. Each 380 x 380 nm 2 map shows 



a pattern of stripes parallel to H . At 0.25 T, six bright 
stripes are seen, with no discernible alternation in ampli- 
tude. At higher H, the stripes increase in number and 
are largely unchanged in width, while amplitude alterna- 
tion becomes noticeable. As H is increased to 0.6 T, the 
major stripes dominate over the minor stripes. The mi- 
nor stripes fade further at 0.8 T and are not discernible 
above 1 T, effectively doubling the stripe periodicity. 

To quantify the stripe periodicities and amplitudes, we 
Fourier transform the Gn(r) data and plot the normalized 
Fourier intensity |Gn(fe)| in Fig. 2(b). In the insets, the 
|Go(fc)| maps show prominent Fourier peaks evolving in 
both fc-space location and intensity. In the main panels, 
the intensities of these peaks are shown by line traces 
of |Go(fc)| from the origin through the Fourier peaks. 
Looking across the panels in Fig. 2, we associate the 
two prominent Fourier peaks with the nearest-neighbor 
(NN) and next-nearest-neighbor (NNN) stripe periodici- 
ties. The NN peaks arc approximately twice as far from 
the fc-space origin as the NNN peaks, indicating that the 
major and minor stripes in Go(r) appear at regularly- 
spaced intervals. At low fields, the NNN peak generally 
has lower intensity than the NN peak. However, as the 
field reaches 0.6 T and above, the relative peak intensities 
are reversed, consistent with fading of the minor stripes 
and doubling of the dominant stripe periodicity. 

To visualize the reversal of NN and NNN peak inten- 
sities, we analyze how the Fourier peaks evolve over the 
entire field range of our measurement. Fig. 3(a) shows a 
color plot of the line traces of |Go(fe)|, for fields between 
0.08 T and 6.5 T applied along [100]. It can be seen that 
the NN peak dominates at low fields and the NNN peak 
at high fields, with a shift in their relative intensities oc- 
curring near 0.7 T. The k positions of both peaks appear 
to scale with VH, as indicated by the dotted curves. We 
note that the intensity shift occurs only for H || [100] but 
not H j| [110]. For the latter field orientation, the NN 
peak dominates the NNN peak both below and above ~ 
0.7 T. This dependence on in-plane field direction is fur- 
ther discussed below, in the context of thermal transport 
data and band-structure and energy-gap anisotropics. 

From the field evolution of the Fourier peaks, we obtain 
information about the in-plane VL. Fig. 3(b) shows a log- 
log plot of the stripe periodicity versus field, which is the 
main result of this Letter. The data points correspond to 
the dominant Fourier peak at each field, for H || [100]. It 
is clear that the data points fall along either of two lines, 
and show a discrete shift near 0.7 T, which we denote 
as H* . This periodicity shift is sharp, namely ±0.1 T 
as defined by the reversal of the relative intensities of 
the Fourier peaks, and shows no intermediate values of 
the stripe periodicity. By fitting the data points to the 
function /3H a , we find that a is -0.51(5) for H < H* 
and -0.46(5) for H > H* , both in good agreement with 
the expected 1 / VH dependence of the VL spacing. On 
the other hand, we find that (3 (in units of nm T -a ) is 
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FIG. 2. (color online), (a) Spatial maps of the normalized zero-bias conductance Go(r), taken on the afe-surface of 2_ff-NbSe2 
at 300 mK with field applied along [100] (black arrow). Each 380 x 380 nm 2 map shows a pattern of stripes parallel to the field. 
The data for 0.25 T is reproduced from Ref. 42. (b) Fourier intensity |Go(fe)| computed from each map, with the inset showing 
|Go(fc)| in the positive quadrant of fc-space, and the main plot showing a line trace of Go(fc)| from the origin through the Fourier 
peaks. To facilitate comparison across fields, the intensity is normalized by the highest peak in each plot. The alternating stripe 
amplitude gives rise to two Fourier peaks corresponding to the nearest-neighbor (NN, solid arrows) and next-nearest-neighbor 
(NNN, open arrows) stripe periodicities. As the field is increased to 0.6 T and higher, the relative intensities of the NN and 
NNN peaks are reversed, consistent with fading of the minor stripes that effectively doubles the dominant stripe periodicity. 



respectively 42(3) and 92(6), indicating a scale change 
in the VL parameters. It is convenient to express these 
results in terms of the Abrikosov VL parameter ao = 
(4/3) 1 / 4 -\/0o7s for a hexagonal lattice, where 4>o is the 
flux quantum and B is the magnetic induction. For H \\ 
[100], the stripe periodicity is thus 0.86(7)ao below H* 
and 1.9(l)ao above H*. 

To interpret these stripe periodicities terms of a dis- 
torted VL, we consider the theory of anisotropic super- 
conductors in the mixed state. According to Ref. 43, an 
in-plane field would produce a distorted hexagonal VL of 
either the A-type or B-type shown in Fig. 4, with the 
respective basis vectors: 

b[ A) = a c/V7, b[ A) = a (6/V7 + V^70)/2 

b[ B) = aoVlV, b[ B) = a a {^ihc + ^yy)/2 (1) 

where y is parallel to the surface and perpendicular to 
the field, and 7 is an anisotropy parameter [44]. Sur- 
face projections of these vectors can be compared to the 
stripe periodicities we observed [45] . Prior studies of 2H- 
NbSe2 have only seen the B-type VL: Gammel et al. [46] 
observed 7 = 3.2(1) using neutron diffraction at 5.2 K 
and 0.8 T for H || [110]; Hess et al. [39] observed 7 of 
up to 3.3 using STS at 300 mK and off c-axis fields up 
to 0.5 T. Our STS data for H || [100] below H* and 
H || [110] are in good agreement with these prior results. 
That is, by associating our NN stripe periodicity with 
surface projections of the first two rows of vortices for 



the B-type VL, we get 0.86(7)a = b\ ' ■ y = Vf } ■ y/2 
with 7 = 3.0(5). For H | [100] above H* , a regime for 
which there have been no prior measurements, our data 
cannot be easily reconciled with the B-type VL. That is, 
in order to stretch b\ ' ■ y from 0.86(7)ao to 1.9(l)ao, 7 
has to increase from w 3 to 14, which is much larger than 
expected (see Supplemental Material [36] for a discussion 
of 7 in 2#-NbSe 2 ). 

To explain the large increase in stripe periodicity for 
H || [100] above if*, we propose a reorientation transi- 
tion of the B-type VL to the A-type VL. In this scenario, 
a similar geometric analysis for the A-type VL, namely 
letting 1.9(l)ao = b% yields 7 = 4.8(5), which is more 
consistent with the available data on 2ff-NbSe2 [39, 46]. 
For the A-type VL, it is important to note that the vor- 

(A) 

tices along b\ are collinear with the tip and thus do not 
show up twice on the surface projection. This geomet- 
ric difference effectively doubles the periodicity between 
adjacent stripes as compared with the B-type VL. Thus, 
the reorientation transition between these two types of 
distorted hexagonal VL provides a natural explanation 
for the stripe periodicity shift. It should also be em- 
phasize that the in-plane VL reorientation transition we 
observed is first-order in character, as evidenced by the 
discreteness of the periodicity shift. This observation is 
in contrast to the second-order VL transition observed 
in MgB2 [16], which involves the gradual rotation of an 
undistorted hexagonal VL for if || c-axis. 

Our results can be directly correlated with the multi- 
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FIG. 3. (color online). Fourier analysis of the field evolution 
of Go(r). (a) Line traces of |Go(fc)| for fields between 0.08 
T and 6.5 T applied along [100]. The two dashed curves 
are a guide to the eye, indicating square-root scaling of the 
peak positions. The reversal of NN and NNN peak intensities 
near H*= 0.7 T is associated with fading of the minor stripes 
and doubling of the dominant stripe periodicity, (b) Log-log 
plot of the stripe periodicity corresponding to the dominant 
Fourier peak at each field. The two dashed lines are power- 
law fits to the data below and above H* , and the dotted line is 
the lattice parameter for a hexagonal Abrikosov lattice. The 
sharp shift near H* shows no intermediate values of the stripe 
periodicity, consistent with a first-order transition. 



gap characteristics seen in heat capacity [28] and thermal 
transport [23] measurements on 2i/-NbSe2. First, the H* 
we observed coincides with the in-plane field at which 
both the Sommcrfcld coefficient and thermal conductiv- 
ity change slope. These characteristics have been inter- 
preted as a distinct signature of multigap pairing, similar 
to the case of MgB2 [20] , whereby field-driven collapse of 
a weak gap as the field is increased above H* causes a de- 
crease in the overall depairing rate. Second, a hysteretic 
anomaly is seen in the thermal conductivity at 380 mK, 
occurring near H * for current perpendicular but not par- 
allel to the in-plane field [23]. The hysteretic behavior 
is consistent with a first-order transition, and although 
the in-plane field direction was not specified in Ref. 23, 
the directional dependence of this anomaly is remarkably 




FIG. 4. Two types of distorted hexagonal vortex lattices for 
anisotropic superconductors with the magnetic field applied 
in-plane. The B-type and A-type lattices are sketched for the 
anisotropy parameter 7 being 3.0 and 4.8, respectively. Our 
STS technique essentially measures the lengths of the basis 
vectors bi and f>2 projected onto the sample surface. Our data 
for H || [100] is consistent with a reorientation transition from 
B-type to A-type, which effectively doubles the predominant 
periodicity in the observed stripe pattern. 



similar to what was also seen in our data. Thus the VL 
reorientation we observed provides a natural explanation 
for the unusual phase transition reported by Ref. 23. 
These experimental correlations indicate a VL reorienta- 
tion transition that is driven by field depairing over parts 
of a multi-sheeted Fermi surface having both uniaxial and 
in-plane anisotropics. This scenario is plausible consid- 
ering the multi-sheeted and multi-gapped Fermi surface 
known to exist in 2/f-NbSe 2 [26, 27, 33, 34]. Detailed 
theoretical study is needed to understand how the ac- 
tual band structure and superconducting-gap anisotropy 
of 2i?-NbSe2 affect its in-plane VL, following a similar 
analysis for MgB2 [47]. 

In the context of multiband pairing, we can also con- 
sider an alternative model for the stripe periodicity shift 
that involves a change in the vortex current distribution. 
Essentially, by regarding the vortex current as consisting 
of multiple superfluid components [48], we can relate its 
radial distribution to the dimensionalities of the Fermi 
surface sheets that take part in the pairing. For 2H- 
NbSe2 above H* , field suppression of pairing on the 3D 
sheet [22, 28] would make the vortices less extended along 
c and limit our ability to detect the second vortex row, 
thus doubling the dominant stripe periodicity. The in- 
plane VL could then remain in the B-type orientation 
above H* , without altering 7, and still produce the peri- 
odicity shift observed. However, to validate this model, 
it would be necessary to calculate the vortex current dis- 
tribution taking into account the actual band structure 
of 2i7-NbSe2 . It would also be necessary to explain why 
the periodicity shift we observed occurs for H || [100] but 
not [110]. 

In summary, using cryomagnetic STS to laterally im- 
age subsurface vortices, we have observed a first-order 
transition of the in-plane VL in 2i7-NbSe2 near 0.7 T 
at 300 mK. The transition is characterized by a discrete 
shift in the projected spacing of a distorted hexagonal 
VL, consistent with its field-driven reorientation. The 
field at which the transition occurs and its dependence 
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on in-planc field direction are correlated with the multi- 
gap characteristics seen in both heat capacity and ther- 
mal transport. These results indicate that the multiband 
pairing of 2iJ-NbSe2 plays an important role in the VL 
transition we observed. 
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FIG. SI. Atomically-resolved topographic image (7.5 x 7.5 nm 2 ) of 2f/-NbSe2 at 300 mK showing 
hexagonally arranged Se atoms and charge-density-wave corrugation, both also visible in the Fourier 
transformed image (inset). 



S2. Anisotropy parameter of 2_£/-NbSe 2 



In the context of a superconductor with uniaxial anisotropy, the anisotropy parameter 7 
is defined as the ratio of superconducting length scales along the c-axis and in the afe-plane. 
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Early works considered 7 to be related to the effective mass anisotropy [1]. In the case of 
multiband superconductors such as MgB 2 , recent experiments [2-4] and theory [5] suggest 
that 7 is in general a function of the temperature T and magnetic field amplitude H . Since 
2i/-NbSe 2 has considerable Fermi velocity (vp) and superconducting gap (A) anisotropy in 
the afr-plane, its 7 can also depend on the in-plane field orientation. Here, we review previous 
measurements on 2iJ-NbSe2 and discuss the expected 7 for 300 mK and H || [100] above 
H*, a regime that has not been previously studied. 

2if-NbSe2 is a superconductor with 3D, quasi-2D and 2D Fermi surface sheets [6]. At 
T <C T c , all the bands are involved in the pairing [7] and thus determine 7. At H > H*, 
however, 7 is expected to increase as pairing on the 3D sheet is suppressed, a phenomenon 
also seen in MgB 2 [8]. 7 can be approximated as J\(T) = A|j c /Aj_ c near H c i and ■Jh c2 (T) = 
H C 2,±c/H C 2 t \\ c near H c2 , where A is the London penetration depth, if cl and H c2 are the lower 
and upper critical fields respectively, and the direction is denoted by the subscripts [5]. 

Previously, Hess et al. [9] measured 7 of up to 3.3, using scanning tunneling spectroscopy 
(STS) at 300 mK and off c-axis fields of up to 0.5 T, below H* = 0.7 T. Above H*, due 
to field suppression of pairing on the 3D band [7, 10], we expect a higher 7 than the one 
measured by Hess et al. [9]. Gammel et al. [11] measured 7 = 3.2(1) using neutron diffraction 
at 5.2 K and 0.8 T with H || [110]. However, since T ~ T c in the latter experiment, thermal 
mixing renders the system less anisotropic. Using the temperature dependence of ^h c2 [12] 
we extrapolate a ~20% increase in 7 at 300 mK from the value of Gammel et al. [11]. 

We now discuss the dependence of 7 on in-plane field orientation. At H = 0.5 T, Hess 
et al. [9] found a 7 that was ~ 10% higher for H || [100] than for H || [110]. In light of recent 
measurements of the Fermi surface of 2if-NbSe2 [13], this can be understood as coming from 
the sixfold anisotropy of the 2D and quasi-2D bands in the afe-plane, which have vp minima 
and A maxima along Y — K, or [110] in real space. Qualitatively, this yields £[100] > £[110], 
where £ = hvp/nA is the coherence length. Thus, i? C 2,[ioo] > #c2,[no] , where i? C 2,[ioo] = 
0o/27r£[ooi]£[iio], #c2,[ho] = 0o/27r£[ooi]£[ioo] and O is the flux quantum. The anisotropy ^ Hc2 
is then larger for H || [100] than for H || [110]. This effect is particularly significant at 
H > H* where the 2D and quasi-2D bands dominate. 

The results of the above analysis are consistent with the field dependence of 7 deduced 
from our data. For H \\ [100] above H*, we observed an A-type vortex lattice with 7 = 
4.8(5). For H \\ [110] both below and above H*, we observed a B-type vortex lattice with a 
lower 7 of 3.0(5). It is worth noting that a similar analysis has been performed for both 
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and 7# c2 in MgB 2 [14], but to our knowledge no such work has yet been carried out for 2H- 
NbSe2- To further elucidate the effect of in-plane field orientation on the superconducting 
anisotropy of 2f/-NbSe2, it would be vital to reexamine previous measurements of its H C 2 
[12], Sommerfeld coefficient [10] and thermal conductivity [15] under known in-plane field 
orientation. 
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